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KiloHertz quasi— periodic oscillations difference frequency 
exceeds inferred spin frequency in 4U 1636—53 
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ABSTRACT 

Recent observations of the low-mass X-ray binary 4U 1636-53 with the Rossi X- 
ray Timing Explorer show, for the first time, a kiloHertz quasi-periodic oscillation 
(kHz QPO) peak separation that exceeds the neutron star spin frequency as inferred 
from burst oscillations. This strongly challenges the sonic-point beat frequency model 
for the kHz QPOs found in low-mass X-ray binaries. We detect two simultaneous 
kHz QPOs with a frequency separation of 323.3±4.3 Hz in an average Fourier power 
spectrum of observations obtained in September 2001 and January 2002. The lower 
kHz QPO frequency varied between 644 Hz and 769 Hz. In previous observations of this 
source the peak separation frequency was —250 Hz when the lower kHz QPO frequency 
was ~900 Hz. Burst oscillations occur in 4U 1636-53 at —581 Hz and possibly at half 
that frequency (290.5 Hz). This is the first source where the peak separation frequency 
is observed to change from less than (half) the burst oscillation frequency to more than 
that. This observation contradicts all previously formulated implementations of the 
sonic-point beat frequency model except those where the disk in 4U 1636-53 switches 
from prograde to retrograde. 
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1 INTRODUCTION 

In recent years kiloHertz quasi-periodic oscillations (kHz 
QPOs) have been found in power spectra of more than 
20 low-mass X-ray binaries (LMXBs) using observations 
mad e with the Rossi X -ray Timing Explorer (RXTE) satel- 
lite (van der Klis 2000). In most sources the kHz QPOs are 



low frequencies were foun d with the EXOSA T and GINGA 
satellites (for a review see van der Klis 1989| ). Observations 
with the RXTE satellite have also reve aled low-frequency 
QPOs in a number of ato ll sources (e.g. 
Ford & van der Klis 19981) 



Homan et al. 199£ 



found in pairs. The high frequencies of these QPOs suggest 
that at least one of the two peaks reflects the orbital motion 
of matter close to the neutron star. Another high frequency 
pheno menon was detected with RXTE; during some type I 



The atoll source 4U 1636-53 is one of the few sources 
that displays all of the high-frequency QPO phenomena that 
have been observed in LMXBs; two kHz QPOs ha ve been 



observed (Zhang et al. 1996; Wijnands et al. 1997) as well 



as a sideband to the lower o f the two kHz QPOs (Jonker 



X-ray Lists nearly mhpw. oscillations wWh .lightly drift Mendez fc van der Klis 2000|) and burst oscil ations at -581 



in frequency, the so-called burst oscillations, have been ob- 
served in the pow er spectra of ten sources (for a review see 
Strohmayer 2001). The X-ray flux-modulation is consistent 
with being due to the changing aspect of a drifting hot spot 
on the surface of the neutron star. Therefore, these burst 
oscillations are thought to reflect the spin frequency of the 
neutron star. In high luminosity 'Z'-source LMXBs QPOs at 



Hz (Zhang et al. 1997; Strohmayer et al. 1998). Miller (1999) 



presented evidence that these oscillations are in fact the 
second harmonic of the neutron star spin of —290.5 Hz. 
Howeve r, using another d ataset these findings were not con- 
firmed (Strohmayer 2001). 

Since the frequency difference between the two kHz 
QPO peaks (the peak separation, Av) is nearly equal to 
(half) the burst oscillation frequency a be at frequency mech- 
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anism was proposed for the kHz QPOs (Strohmayer et al 
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1996). In such a model the higher frequency QPO is at- 
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tributed to the orbital frequency of clumps of plasma at 
a special radius near the neutron star while the lower fre- 
quency QPO is due to the beat between the orbital frequency 
and the stellar spin frequency. A specific model incorporat- 
ing the beat frequency mechanism, the sonic-point beat fre- 
quency model, is now one of the leading models exp laining 
kHz QPOs in LMXBs ( |Miller, Lamb fc Psaltis 1998| ). How- 
ever, subsequent findings complicated the picture. As was 
first shown in Sco X-l, the peak separation decreases as th e 
kHz QPO frequencies increase ( van der Klis et al. 1997 ). 
Furthermore, in 4U 1636-53 the peak separation w as found 



to be le ss than half the burst oscillat ion frequency ( Mendez 



van de: Klis & van Paradijs 1998). These findings could 
be explained within the sonic point beat frequency model 



and a time resolution of 128 s to visualize the time evolu- 
tion of the lower kHz QPO for all observations except for 
the power spectra of observation 60032-01-21-00 which were 
rebinned to a time resolution of 162s (segment 7) and 256s 
(segment 8), and the power spectra of observation 60032- 
01-22-00, which were treated separately. We traced the fre- 
quency of the lower kHz QPO and fitted the power spectrum 
in a window of 200 Hz centered on the traced frequency with 
a Lorentzian and a constant. All peaks in each power spec- 
trum of 128 s, 162 s, and 256 s had a significance of ~ 3cr or 
more. The frequency of the bin associated with the peak of 
the fitted Lorentzian was taken to be the frequency of the 
QPO in that power density spectru m. Finally, we used the 
shift-and-add meth od described by Mendez, van der Klis & 



( Lamb fe Miller 200l|) by taking into account the effect of van Paradijs (1998) to shift the lower kHz QPO to the same 



the inward velocity of the plasma on the emerging kHz QPO 
frequencies. A firm prediction of this revised sonic-point 
beat frequency model is that for prograde spinning accre- 
tion disks the kHz QPO peak separation is smaller than the 
neutron star spin frequency. This follows from the detailed 
description of the plasma flow patterns in this model which 
produce a number of corrections to the observed frequencies 
relative to the true orbital and be at frequency, all of w hich 
are expected to make Av smaller (Lamb & Miller 2001). In- 



deed, until now, in various sources the peak separation was 
always found to be consistent with, or less than (half) the 
burst oscillation frequency. 

In this Letter we show that in recent observations of 
4U 1636-53 the kHz QPO peak separation is significantly 
larger than h alf the burst oscillation frequency. Using other 
observations Mendez, van der Klis & van Paradijs (1998) 



had shown earlier that the peak separation was significantly 
smaller than half the burst oscillation frequency. This means 
that in order for the sonic-point beat frequency to explain 
both observations further changes will have to be made in 
the description of the flow pattern. The only change in the 
flow pattern previously described that could produce the 
observed change in the separation frequency is that in which 
the accretion disk changes from prograde to retrograde when 



the lowir kHz QPO moves from ~750 Hz to ~800 Hz. 



2 OBSERVATIONS AND ANALYSIS 

4U 1636-53 was observed in September 2001 and January 
2002 with the proportional counter arr ay (PCA; lahoda 
et al. 1996|) onb oard the RXTE satellite ( |Bradt, Rothschild 
fc Swank 1993). Due to the reduced duty cycle of the PCA 



detectors only a subset of the 5 detectors is operational for 
most of the time. In this analysis we only used data in which 
4 of the 5 detectors were operational. A log of the observa- 
tions is given in Table [l} In total ~28.4 ksec of data were 
used. Data were obtained using an Event-mode in which in- 
formation on detected X-ray photons is stored onboard and 
telemetered to the ground on an event-to-event basis. The 
time resolution of the data is ~125/^s and RXTEs effective 
energy range (2-60 keV) is covered by 64 channels. 

Using this data we calculated power density spectra of 
segmen ts of 64 s up to a Nyquist frequency of 4096 Hz in 



(arbitrary) frequency in each power spectrum before aver- 
aging. This assures that at least for timescales longer than 
256 s the determined peak separation is not influenced by 
the fact that the tw o kHz QPOs have a different freque ncy- 
amplitude relation (Mendez, van der Klis & Ford 2001). 

We averaged the aligned power spectra and fitted the 
average power spectrum from 256-1500 Hz with a constant, 
to represent the power introduced by the Poisson counting 
noise, and two Lorentzians to represent the two kHz QPOs. 
Errors on the fit parameters were determined using A^ 2 = 
1.0 (la single parameter). The observation 60032-01-22-00 
was treated separately since we were not able to find a more 
than 3<r lower kHz QPO using data stretches shorter than 
or equal to 256 s. Instead we averaged the power spectra 
in this observation without shifting the data. We note that 
by doing so we may have artificially introduced deviations in 
the kHz QPO peak separation if the two kHz QPOs moved in 
frequency during the observation for reasons outlined above. 

The fits of the average power spectra were good, with a 
reduced x 2 slightly larger than 1 for 146 degrees of freedom. 
We found that, combining all observations, except 60032-01- 
22-00, the frequency separation between the two kHz QPO 
peaks was 323.3±4.3 Hz, significantly larger than 29 0.5 Hz 
the put ative spin frequency according to an analysis of 



Miller 



(1999), but still significantly smaller than 581 Hz, the fre- 



quency of t he strong burst oscillation s found by ^hang et al 
(1997) and [Strohmayer et al. (1998)| (see Figure ||). We then 
divided the data set into three parts according to the fre- 
quency of the lower kHz QPO and measured the peak sep- 
aration in each of the parts. For each part we also averaged 
the power spectra without shifting in order to measure the 
average frequency of the lower kHz QPO and that of a low- 
frequency QPO changing in frequency from ~24-42 Hz. In 
Table |^ we present our results on the kHz QPOs. The frac- 
tional rms amplitude of the low-frequency QPO changed 
from 7.8±0.6 per cent to 4.3±0.7 per cent (2-60 keV) while 
its frequency increased from 23.6±2.4 Hz to 42±4 Hz and 
its FWHM was consistent with being constant at 26±9 Hz. 
In Figure |^ we plot the measured peak separation (Au) as 
a function of the frequency of the lower and upper kHz 
QPO (left and right panel, respectively). The squares are 
our new measurements using the shift-and-add method, the 
diamond is our measurement using observ ation 60032-01-22- 
00, and the dots are the measurements of Di Salvo, Mendez 



one energy band covering the total PCA energy range. We 



van der Klis (200g_ 



used a dynamical power spectrum displaying consecutive 1996-2000 (see also IWijnands et al. 1997 ; Mendez, van der 



2J ba; 
lWijn. 



based on RXTE data obtained from 



power spectra rebinned to a frequency resolution of 2 Hz Klis fc van Paradijs 199£). Clearly, Av changes from more 
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Table 1. Log of the observations used in this analysis. 



Observation ID 


Used segments 


Amount of good data (ksec) 


Count rate" cnt/s/PCU 


60032-01-09-00 


1 


3.3 


230 


60032-01-09-03 


all 


2.5 


220 


60032-01-19-00 


1, 2, 3 


8.6 


205 


60032-01-21-00 


7, 8 


6.7 


180 


60032-01-22-00 


1, 2 


7.3 


180 



a Average background subtracted count rate per PCU (2-60 keV) 



than h; .If the burst oscillation frequency to less than that. 



The average Av of the three leftmost points is 327.2±4.2 
Hz and that of the dots is 247.1±1.9 Hz making the jump 



in Av 73.1 ±4.6 Hz. 



der Klis fc Jonker 2001 ) could also bias Av measurements. 
However, those ~5-6 Hz QPOs are only found when atoll 
and Z sources are in a high-luminosity state (e.g. Middled 



itch fc Priedhorsky 1986; Wijnands, van der Klis & Rijkhorst 



3 DISCUSSION 

We found that in 4U 1636-53 the kHz QPO frequency sep- 
aration (Av) changes as a function of the frequency of the 
lower kHz QPO from 323.3±4.3 Hz to 242.4±3.6 Hz (Fig. | 
left); the maximum peak separation is significantly larger 
than 290.5 Hz (half the burst oscillation frequency of ~581 
Hz; see Section |l|), while the minimum peak separation is 
significantly smaller than that. The change in Av as a func- 
tion of upper kHz QPO frequency is abrupt; the data are 
consistent with a sudden jump at v — 1050 Hz (Fig. |^ right). 
We note that similar jumps in Av as a function of the upper 
kHz QPO have also been seen in 4U 1728-34 and 4U 1608- 
52, although in those cases Av always remained smaller or 
equal to (half) the burst oscillation f requency (Mendcz fc 
van der Klis 19991; |Mendez et al. 1998| ; [Lamb fc Miller 200l[ 
Muno et al. 200l|). This is the first time that the kHz peak 



separation has been shown to be significantly larger than 
the inferred neutron star spin frequency and also the first 
time that Av has been seen to vary between less and more 
than (half) the burst oscillation frequency. 

In previous work burst oscillations detected in the dip- 
ping burster 4U 1916-05 at ~272 Hz also seemed to occur at 



too low a frequency for Av (Galloway et al. 2001). However, 
in that source the kHz QPOs are weak and sometimes very 
broad, therefore whole segments of several ksec of dat a had 
to be averaged to obtain a detection ( Boirin et al. 200o| ) . This 
means that there is no certainty that the kHz QPOs were 
detected simultaneously or that their apparent separation 
was not affected by frequency shifts combined with changes 
in their power ratio. Therefore, the measured Av may not 
reflect the true peak separation. In the current work, using 
the shift-and-add method we excluded the possibility that 
effects changing Av on timescales longer than 256 s have in- 
fluenced the observed Av. On timescales shorter than 256 s, 
changes in the rms amplitudes could bias the peak sepa- 
ration for reasons outlined above. However, since our mea- 



suremc its show that the lower kHz QPO rms amplitude 
increased while the rms amplitude of the upper kHz QPO 
decreased this effect will have made the observed Av smaller 
than the true peak separation assuming the rms amplitude- 
frequency relation is the same for timescales shorter than 
256 s. Fast changes in the upper kHz QPO frequen cy, asso- 
ciated with the 5-6 Hz QPO in Sco X-l (~20 Hz; |Yu, van 



199£). Changes in the kHz QPO frequencies associated with 
milliHz QPOs present in 4 U 1636-53 and 4U 160 8-52 in 
a small luminosity inte rval (Revnivtsev et al. 2001 ) are so 
small (less than a Hz; Yu & van der Klis 2002) that it is 
unlikely that they alter our findings significantly, although 
we note that such a milliHz QPO was present during some 
of the observations. 

Systematic variations in the profile of the kHz QPOs 
could in principle cause the kHz QPO separation frequency 
to become large at low kHz QPO frequencies and small at 
high frequencies since the centroid frequencies have been 
measured assuming a Lorentzian peak profile. Since our fits 
were good (reduced \ 2 ~1) an d no systematic residuals were 
apparent, we estimate that shifts in the centroid frequency 
due to asymmetries in the peak profiles could at most be 
as large as the FWMH of the peak divided by its signifi- 
cance, i.e. negligible in case of the lower kHz QPO but up to 
~10 Hz in case of the upper kHz QPO. For power spectra 
with a lower kHz QPO in the frequency range 650-750 Hz, 
the average data— model residuals within one FWHM are 
less then 0.15 per cent of the power in the peak. Further- 
more, there is no evidence for a sudden change of the QPO 
profile which could produce the observed abrupt decrease in 
Av. When using Gaussians to model the kHz QPO peaks the 
fits were also good (reduced \ 2 ~1) an d we obtained results 
consistent with those using Lorentzians. We conclude that 
it is unlikely that our results can be explained by changing 
asymmetries in the peak profile of the kHz QPOs. 

The principal motivation for a beat frequency model is 
the relative closeness of the inferred spin frequency to Av. 
We found that Av is consistent with being distributed sym- 
metrically around half the burst oscillation frequency, em- 
phasizing the apparent connection between the inferred spin 
and Av. The present formulation of the sonic-point beat- 
frequency model can explain the obs erved decrease in Av a s 
a function of kHz QPO frequency flLamb fc Miller 200l[ ). 
However, all versions of the sonic-point beat-frequency 
model described so far involving a prograde spinning ac- 



cretion disk (Miller, Lamb & Psaltis 199£; Lamb & Miller 



2001) predict that Av is always less than or equal to the spin 
frequency of the neutron star for Keplerian orbital frequen- 
cies larger than the neutron star spin frequencies; since our 
results for 4U 1636-53 show that Av can be larger than the 
neutron star spin frequency this prediction is now falsified. 

Strohmayer (2001) showed that it is still unclear 
whether the burst oscillations at ~581 Hz reflect the spin 
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Table 2. The kHz QPO frequency separation (Ai/), the frequency, the fractional rms amplitude (2-60 keV), and the full-width-at-half- 
maximum (FWHM) of the lower kHz QPO (ui ower kHz)> the fractional rms amplitude (2-60 keV) and the FWHM of the upper kHz 
QPO {u upper kHz)- 





Lower kHz QPO 






Upper kHz QPO 


Peak separation 


rms (%) 


FWHM (Hz) 


"lower kHz (Hz) a 


rms (%) 


FWHM (Hz) 


At/ (Hz) 


7.7±0.4 


52±8 


644.2±3.2 


6.9±0.5 


52±11 


326.9±5.3 


7.6±0.3 


15.7±2.3 


687.8±2.1 


6.8±0.9 


74±35 


325.1±11.5 


8.5±0.1 


8.1±0.3 


723.0±0.9 


5.6±0.7 


83±31 


329.5±8.9 


8.7±0.1 


7.2±0.7 


769.2±1.2 


5.1±0.9 


139±62 


288.6±31 



a The frequency of the lower kHz QPO (ui ower kHz) an d an the parameters on the first row were determined without shifting the data. 
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Figure 1. Average power density spectrum combining all shifted data used in this paper. The peak of the lower kHz QPO is off the 
vertical scale. The kHz QPO frequency separation is 323.3±4.3 Hz. The power spectra have been shifted before averaging, therefore, only 
the frequency difference between the peaks is meaningful. 
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Figure 2. Left: The kHz QPO frequency separation as a function of the frequency of the lower kHz QPO. Right: The kHz QPO 
frequency separation as a function of the frequency of the upper kHz QPO. The frequency of half the burst oscillation frequency (290.5 
Hz) is indicated by the dashed line. The peak separation varies from significantly more than half the burst oscillation frequency [diamond 
(unshiftcd data) and squares (shifted data); this work] to significantly less than that (dots; Di Salvo et al. 2002; see also Wijnands et al. 
1997 and Mendez et al. 1998). Error bars on the frequency of the lower kHz QPO are generally smaller than the size of the symbols. 
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frequency of the neutro n star or its second harmonic as sug- 
gested by Miller (1999). Therefore, one might argue that the 



spin frequency is really at 581 Hz (not 290.5 Hz) and hence 
that Av is less than the spin frequency. This would, how- 
ever, take away the principal motivation of beat -frequency 
models, since the neutron star spin frequency and the kHz 
QPO peak separation are in that case completely different. 

That in the sonic-point beat-frequency model the neu- 
tron star spin and Av are not exactly equal is explained by 
the way in which the beat-frequency interaction produces 
the observed frequencies. The key element is the r adial mo- 
tion of the plasma clumps near the sonic radius (Lamb & 
Miller fo0l| ). This makes the observed lower kHz QPO fre- 
quency somewhat larger than the beat-frequency since it 
squeezes the spatial separation between peaks of enhanced 
mass flow. It also makes the upper kHz QPO frequency 
smaller than the orbital frequency since it makes the clump's 
footpoint move backward in a corotating frame. It might be 
possible to find other eff ects in the plasma flow patterns be- 
sides those discussed by Lamb & Miller (2001) (all of which 



make Av smaller) which may lead to larger values of Av at 
lower kHz QPO frequencies. Fully relativistic simul ations of 
the gas flow and radiation transport reported by 



Lamb & 



Miller I 2001) have not, so far, shown evidence for this. If, 
however, the sense of rotation of (part of) the accretion disk 
would change from prograde to retrograde when Av changes 
from less to more than the spin frequency, then our results 
could still be in accordance with current formulations of the 
sonic-point beat-frequency model. This also requires that in 
our observations for Av larger than 290.5 Hz the lower kHz 
QPO reflects the (near) Keplerian orbital frequency, whereas 
for those points in Figure H where Av is smaller than half 
the burst oscillation frequency, the Keplerian orbital motion 
is represented by the upper kHz QPO. 

We detected a low-frequency QPO between ~25-40 Hz 
in 4U 1636-53 simultaneously with the kHz QPO pair. The 
relation between this low-frequency Q PO and the lower kHz 



QPO is similar to the one found by Ford fc van der Kli , 
(1998)|]for 4U 1728-34. It has been suggested that such 
low-frequency QPOs in atoll sources are related t o the Z 
source Hor i zontal Branch Oscillations (HBO; e. g. Homai 
et al. Ip9§; |Psaltis, Belloni fc van der Klis 1999|). Accord- 



ing to the magnetospheric beat -frequency model these HBO 



origi nate near the mag netospheric radius (Alpar & Shaham 



1985; L imb et al. 1985). The frequency of the low-frequency 
QPO in 4U 1636-53 increased as the frequency of the lower 
kHz QPO increased. Therefore, if the disk near the sonic 
radius is indeed retrograde within our observations the ret- 
rograde annulus must extend at least up to the magneto- 
spheric radius if the sonic-point and magnetospheric beat- 
frequency model both apply. Partially retrograde spinning 
disks have been proposed as an ex planation for observed 
torque reversals (Bildsten et al. 1997) in slowly rotating neu- 



tron st ars with a dipole magneti c field strength of 10 11 J " 
Gauss (van Kerkwijk et al. 1998). Modelling by van Kerk- 



wijk et al. (1998) shows that the accretion disk can get incli- 
nations of more than 90 degrees. It is not clear if their work 
is also applicable to the case of low-magnetic field LMXBs. 
Clearly, allowing disks with alternating rotation senses com- 
plicates the sonic-point beat frequency interpretation of the 
kHz QPOs considerably. 
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